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We consider the leading post-Newtonian and quantum corrections to the non-relativistic scatter- 
ing amplitude of charged spin-i fermions in the combined theory of general relativity and QED. 
The coupled Dirac-Einstein system is treated as an effective field theory. This allows for a consis- 
tent quantization of the gravitational field. The appropriate vertex rules are extracted from the 
action, and the non-analytic contributions to the f-loop scattering matrix are calculated in the non- 
relativistic limit. The non-analytical parts of the scattering amplitude are known to give the long 
range, low energy, leading quantum corrections, are used to construct the leading post-Newtonian 
and quantum corrections to the two-particle non-relativistic scattering matrix potential for two 
massive fermions with electric charge. 



I. INTRODUCTION 

On the search for a theory of quantum gravity, 
Donoghue il| proposed 12 years ago an interesting new 
way to look at general relativity. He suggested that when 
treating general relativity as an effective field theory Q , 
reliable quantum predictions at the low energies could 
be made, in the same way as chiral perturbation theory 
is used as the low energy approximation of QCD, being 
the effective field theory of QCD. It is well known that a 
field theory need not be strictly renormalizable in order 
to be able to yield quantum predictions at low energies. 
A fundamental quantum theory of gravity does not ap- 
pear in this way, but it is possible to calculate quantum 
corrections order by order in a momentum expansion. 

Having laid the foundations for this new approach, 
Donoghue and collaborators turned their attention to the 
practical applications of this idea. A number of interest- 
ing calculations has been made involving quanturngrav- 
itational corrections to various quantities P, |3j IE Hi ■ 

Prior to the effective field theoretical description of 
general relativity, attempts had been made to find a 
quantum theory for gravity. In particular many proved 
that general relativity was indeed not renormalizable, be 
it pure general relativity or general relativity coupled to 
bosonic or fermionic matter, see e.g. 
course it is a well known fact that general relativity is a 
non-renormalizable theory per se, and these authors suc- 
ceeded in exactly confirming that gravity indeed is explic- 
itly non-renormalizable, with or without matter. How- 
ever, when looked at in the framework of an effective field 
theory, these theories do become order by order renormal- 
izable in the low energy limit. All possible counter-terms, 
also those not present in the initial Lagrangian, are gen- 
erated. When general relativity is treated as an effective 
theory, renormalizability simply fails to be an issue. The 
ultraviolet divergences arising e.g. at the 1-loop level 
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are dealt with by renormalizing the parameters of higher 
derivative terms in the action. Many interesting results 
have been found from this procedure. Most interesting 
for the point of view of this paper is the bosonic quantum 
corrections to the Newtonian/Coulomb potentials (61]. 

When approaching general relativity in this manner, it 
is convenient to use the background field method 0, ■ 
Divergent terms are absorbed away into phenomenolog- 
ical constants which characterize the effective action of 
the theory. The price paid is the introduction of a set of 
never-ending higher order derivative couplings into the 
theory. General relativity thus turns into a minimal the- 
ory in the myriad of higher order terms in the action, 
however still remaining a valid theory for gravitational 
interactions at low energies. The effective action con- 
tains all terms consistent with the underlying symmetries 
of the theory. Perturbatively certain terms of the action 
play leading roles at certain energy scales, hence only a 
finite number of terms are required to be accounted for 
at each loop order. Here we consider the low energy limit 
of the effective theory of quantum gravity. Because we 
work to lowest order, our results are free of the new ad- 
ditional terms that must be appended to the Einstein 
action and which manifest themselves as components of 
the high energy couplings of the effective field theory. 

In ref. the post-Newtonian as well as the quantum 
corrections, that were generated to the Newtonian and 
Coulomb potentials were explicitly found. We wish to 
repeat this calculation, but now in terms of couplings 
to fermions. We wish in particular to see explicitly if 
the post-Newtonian as well as the quantum gravitational 
corrections generated are identical to P| or not. 

We will more or less follow a similar procedure as in 
mostly to avoid confusion about conventions and to make 
it easy to compare the results at the end. However, it is 
by no means a straightforward task to complete a simi- 
lar calculation in terms of fermions. Some obstacles have 
to be overcome compared to bosonic matter in curved 
space, e.g. the issue of introducing fermionic matter into 
curved space-time. Luckily, this issue has been dealt with 
before H E 111 111 IH 111 13 ■ The additional formal- 
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ism required is the introduction of the vierbein formahsm 
and deriving a proper covariant derivative for the spinor 
fields. 

Donoghue devised a particular elegant way to extract 
relevant information in terms of analytical and non- 
analytical contributions to the scattering matrix. This 
was realized through the integrals occurring in the cal- 
culations, and propagation of massless particles. Be- 
cause the post-Newtonian and the quantum corrections 
are fully determined by the non-analytical pieces of the 
1-loop amplitude generated by the lowest order Einstein 
action, it becomes possible to perform this calculation 
completely. We will also only consider 1-loop effects in 
this paper. We will extract the non-analytical parts of 
the full set of 1-loop diagrams needed for the 1-loop scat- 
tering matrix in the combined quantum theory of general 
relativity and QED. As we will see in this paper, and as 
can also seen in 0, 0, H, 0, 13 1 the non-analytical contri- 
butions correspond exactly to the long range corrections 
of the potential. 

We will employ the same conventions as in ref. 0, 
01 . The mostly minus Minkowski metric convention 
(1,-1,— 1,-1) will be used and the natural units are 
{h = c = 1) when nothing else is stated. 

In section [H] we will shortly review the concept of 
vierbein fields and how to introduce a proper covariant 
derivative for the spinors when working with fermions. 
We will also quantize the metric and vierbeins using the 
background field method. We will explain the correspon- 
dence between the metric and vierbein formulation of 
our theory. In section IIIII we will see how to combine 
QED with general relativity by using the vierbein for- 
malism and moreover introduce the ghost fields. Next 
in section IIVI we will focus first on the distinction be- 
tween non-analytical and analytical contributions to the 
scattering matrix amplitude, where after we will define 
the potential. Finally, in the succeeding section we 
will evaluate the Feynman diagrams contributing non- 
analytically to the scattering matrix, in order to con- 
struct the leading corrections to the non-relativistic New- 
tonian and Coulomb potential. We will end this paper 
with a discussion in relation to (sj. In the appendix, the 
vertex rules are presented. 



II. INTRODUCING THE VIERBEIN FIELDS 

At every space-time point xq it is possible to erect a 
set of coordinates ^"^^ locally inertial at the given point 
in question, in accordance with the theory of special rela- 
tivity. This in turn implies that the erected set of coordi- 
nates ^"jj vary from point to point, and that the informa- 
tion about the gravitational field is in fact contained in 
the change of the local inertial coordinate systems from 
point to point. Hence it is possible to express as a 
local function of any non-inertial coordinates (i.e. in a 



general coordinate system) x^ i.e. 



(la) 
(lb) 



evaluating the derivatives at the point of interest, where 
the inverse operation is also given at the point in ques- 
tion. The transformation matrix relating the local in- 
ertial frames to the arbitrary coordinate system is the 
called vierbein field, it is denoted e°'^(x) and is a func- 
tion of x^. 

It is seen that due to the transformations x and 
^ — > a; being nonsingular, ea^{x) is the inverse transfor- 
mation matrix. 

One can find other relations between the vierbein 
fields, e.g 



e fj,eb — Oft 



(2a) 
(2b) 



The metric expressed in general coordinates can be found 
by looking at the proper time in general coordinates 



dr^ ^iMdCix)d^\x) 
where the metric in general coordinates is 
gfiiy{x) = e''^{x)e''i,{x)r]ab 



(3) 



(4) 



where a,b . . . are Lorentz indices and fj,,i> . . . are the gen- 
eral coordinate indices. The inverse of the metric is 



(5) 



g'^-' =ea^{x)eb''{x)r^'''' 
from which a familiar result is obtained 

5^.5"" = (e^e^r;,^,)(ee■^ed"77='') = e%ea^ = 5^ 



Thus the e''^(a;) fields relate the Lorentz axes to the co- 
ordinate axes at each point in space-time, explicitly 



e"^(a; = xq) 



dxi" 



(6) 



When changing the general coordinates from x'^, 
these fields transform as covariant vectors and under lo- 
cal Lorentz transformations they transform as Lorentz 
vectors 



dx" 



(7a) 



^x'^' 

be- (7b) 



The vierbein fields are the only index changing objects 
in this theory. For given covariant/contravariant vector 
fields or also a tensor field, we can refer their components 
at X to the locally inertial coordinate system (x) at xq 
by using the vierbein 



(8) 
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A"' transforms as a collection of four scalars under gen- 
eral coordinate transformations (|7a|) , and under the local 
Lorentz transformations l)7b|l it behaves as a vector. It 
is thus possible to use the vierbeins, to convert general 
tensors into local, Lorentz-transforming tensors, whereby 
shifting the additional space-time dependence into the 
vierbeins. 



A. Gauge transformation of the fermion fields and 
the spin connection 

Using the generators of the Lorentz group aab — 
j[la,lb] for spinors belonging to a given representation 
S{A) = ei^'''''^ab fQj- sniall A in infinitesimal form we get 



for spinor transformations 



Note that the indices commute only with alike (space- 
time-) indices, i.e. Latin characters with Latins and 
Greek characters with Greek. 

This object has exactly the transformation property 
that we were searching for 



, ,ab lab 



ab 



(14) 



where the primed object shows transformation as a ten- 
sor. The covariant derivative for spinor fields now be- 
comes 



(15) 



This will transform as a proper Lorentz spinor and a 
covariant vector under the mentioned gauge transforma- 
tions. 



1 



^' ^{1+ aab\''\x))^ 



(9) 



B. Quantization of the vierbein and metric fields 



Since partial derivatives on fields almost always occur in 
the Lagrangian density functions, it would be worthwhile 
to subject Q to a partial differentiation, and see if it 
transforms covariantly or not. We immediately see from 
the transformation 



1 



1 



^TafcA'"'(x))a^^- 



CJabd^,X'\x) 



^ (10) 



that the fermion field does not transform as a proper 
Lorentz spinor under this operation, even though the 
ordinary derivative is a covariant vector since the 
fermion fields are defined to be scalar objects under coor- 
dinate transformations. In order to get a feasible theory 
(|10|l must transform covariantly, hence we must invent 
a "covariant derivative" for fields that transform in 
this manner, such that in the end result the second term 
in ifTUIl gets canceled. 

We can find this object by studying the behavior of the 
vierbein fields under the transformations H7a(l and (|7bp 
and comparing it with the transformation properties of 
a general four vector . 

The rule for covariant derivative for an object mixed 
in both Lorentz and general coordinate indices A" is 



To quantize the combined theory of QED and gravity 
we will have to quantize the vierbein fields in the same 
way as it is done when quantizing the metric. 

In the background field method the quantum correc- 
tions to general relativity are described by quantum vi- 
brations of the metric tensor, making it possible to ex- 
pand the metric and vierbein fields into two separate con- 
tributions, a classical background field and a quantum 
field 



(16) 
(17) 



where k"^ = 32ttG and the background fields are denoted 
as g^^ and e°^. The quantum part - the graviton field 
- is denoted by /i^^ and c"^, the sum of these being the 
full metric and vierbein respectively. We will soon see 
how the metric and vierbein quantum variables for the 
gravitons are related. Furthermore we find the following 
inverses and other relations 



(18) 



D,Al = d^Al ~ V^^^Al + uj;\Al (11) 

where u)^"'b is called the spin connection. To determine 
the structure of the spin connection one can require that 
the covariant differentiation should commute with the 
operation of index changing, not only just index lower- 
ing/raising , thus using the condition 24] 



(12) 



one finds that the following object does the job 

< = \ (e['"'a[^e^] + e'^''e^'^a[^e,,]e=,,) (13) 



for the vierbeins 



(19) 



for the metric. We have only expanded the inverses to 
first order in the quantum fields, as we need diagrams to 
second order in k. From these relations we see that the 
metric and vierbein quantum fields are related according 
to 



(20) 
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showing us that the quantized metric field is equal the 
quantized symmetric vierbein field to first order in the 
quantum fields, i.e. hfj_^ — Ci^^^ + c^^ = s^^. It is easy to 
deduce the determinants of the vierbein and metric fields 



e = det[e°^] « e(^l 



(21) 



with e — det e""^ , as well as the square- root of the metric 
tensor = yj~det{g^,u) 



-9 



-g 



(22) 



with g = det g^a ■ Finally it will be necessary to expand 
the spin connection in terms of the vierbein fields as well, 
the leading order terms appearing in the calculations are 
(in flat background field) 



Background pi 



(23a) 



w 



First order 
liah 



= 7:dp.aba + 7idbSaf_i - -daSb^i (23b) 



where we have defined a new field, an antisymmetric field 



where m is the mass, Cq '2^ is the electron charge with 
= \eq\ and finally Z?^ = d^ — ieqAf^{x) is the covariant 
derivative. 

To make the action of the Dirac Lagrangian density 
invariant under general coordinate transformations, we 
follow the general procedure, i.e. multiply it with ^/—g 
and at the same time introduce our new covariant deriva- 
tive 



now = 



det(e%) 



e i.e. the determinant of the vierbein 
Finally 7^ = 



to) -0 
(26) 

and we have used 



is the matrix square-root of the metric 

The full generally covariant Lagrangian density includ- 
ing the fermionic degrees of freedom may collectively be 
written as 



--V^g'^^g^'F^^F^f, (27) 



III. THE DIRAC-EINSTEIN SYSTEM AS A 
COMBINED EFFECTIVE FIELD THEORY 

The combined theory of QED in a gravitational field 
is given by the sum of the QED and Einstein Lagrangian 
densities 



This will account for our full theory. The Lagrangian 
density is to be expanded in powers of c^^ (where we 
choose c^jy to be linearly symmetrically equal to h^^, as 
seen earlier) in the case of the Dirac field and only /i^^ 
in the case of the Maxwell fields, specifically we expand 
the Lagrangian density as follows 



avity 



c 



QED 



(24) 



The interacting field theory for Quantum Electrodynam- 
ics is well known, with the Dirac equation minimally cou- 
pled to the electromagnetic field 



£qED — ^Dira 



c 



Maxwell 



{Pitro,, ~ m)^ ~ -^g^'^g^-F^.F.p ^^^^ 



— -^Background H~ ^Linear order 



(28) 



where the ellipses denote second and higher order terms 
that will not contribute at the 1-loop level calculations. 

The Lagrangian density for the photon field can now 
be expanded in powers of /i^^ , explicitly 



^Maxwell = --^hid^^A^d^A^ - d^A^^d^ A^) + 'l^h^'' (d^A.d'^ Ap + d.A^dpA- - daA^dpA" - d^A^d^Ap) (29) 



where the trace of ft- = h°'^ = = h. And likewise for the fermionic part 



(30) 



where the symmetric identity la^^"^ = 5'7a^^'7"^^- theory from the Linear order expansions as stated above 

All the necessary lowest order interaction vertices of in equations and ((201 ■ A summary of these rules is 
fermions, gravitons and photons can be found for the presented in the appendix. 
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In principle we should also expand 

2 

C = e^R (31) 

however it is known that the metric and vierbein for- 
mulations are equivalent for fields with only covariant 
vector indices. The coupling to the vierbein field only 
occur as symmetric combinations of vierbein fields, the 
symmetric combination is as we have seen to linear or- 
der equal to the metric tensor. No new aspects of the 
traditional quantization of the pure gravitational action 
are introduced. We will therefore use the known vertices 
and propagators for the bosonic and gravitational fields. 

We have excluded the antisymmetric fields in all our 
expressions. This is due to the fact that the antisymmet- 
ric fields have propagators that go as ~ k^. This can be 
seen when we fix the gauges in our quantization scheme. 
As pointed out earlier, our theory (i.e. fermions includ- 
ing gravitational effects) has two types of invariances. 
One is the general coordinate transformations l|7a|) . un- 
der which the fermions behave as scalars (since they are 
defined with respect to the local Lorentz frame). The 
other is the local Lorentz transformations (j7b|) . under 
which the fermions transform as spinors. If Einstein ac- 
tion is included, then the coordinate gauge can be fixed 
by choosing the harmonic (de Bonder) gauge 

^'' = -lV^{h,.'''-lK-,,r (32) 

whereas the local Lorentz invariance is broken by choos- 
ing the sum of the squares of the antisymmetric vierbein 
components 

C'^ = -^et,-'a% (33) 

Gauge fixing of both these fields will result in an intro- 
duction of two sets of ghost fields. We do not need to 
be concerned about the ghost introduced due to the an- 
tisymmetric field. In a vierbein description of pure grav- 
ity, the ghosts are never external, furthermore neither 
the antisymmetric vierbein fields nor its ghosts propa- 
gate (they cancel each other 0), thus we will not need 
to calculate vertices for the external ghosts fields. This is 
very reassuring since the pure gravity theory in vierbein 
formulation can be covariantly quantized and is equiva- 
lent to the quantized metric approach. That is we could 
in principle describe the theory without introducing these 
variables. But if we do not have pure gravity and include 
fermions, the antisymmetric fields become coupled to the 
vertices. We need only consider the symmetric fields of 
the interactions. This is due to the fact that we will only 
be interested in the long range corrections to the back- 
ground field, and the antisymmetric fields do not produce 
non-analytic terms to the order at which we are working, 
due to the proportionality factor of its propagator ~ . 
In fact a diagram consisting of at least an antisymmetric 
field and a graviton vertex will at least go as which is 



an order higher than ~ k^. However in a full treatment of 
gravitational interaction between fermionic matter these 
fields will have important contributions. They will most 
likely contribute to higher order calculations. 

IV. THE SCATTERING MATRIX AND THE 
POTENTIAL 

It will be fruitful to make a distinction between non- 
analytical and analytical contributions from the dia- 
grams, in order to compute the leading long range, low 
energy quantum corrections to this theory. 

This distinction originates from the impossibility of ex- 
panding a massless propagator ^ while on the other 
hand we have 



expansion for the massive propagator. No ~ terms are 
generated by the above expansion of the massive propa- 
gator. We see that the non-analytical contributions are 
inherently non-local effects which cannot be expanded in 
a power series in momentum. Thus the non-analytical 
effects appear from the propagation of massless particle 
modes, in our case the gravitons and photons. These non- 
analytic contributions will be governed to leading order 
only by the minimally coupled Lagrangian. 

The non-analytical effects are e.g. terms in the S- 
matrix which go as ~ ln(— g^) or ~ —^=^ the general 

example of an analytical effect is a power series in the 
momentum q. The analytical contributions in the dia- 
grams are local effects, they will always be expandable 
in power series solutions. We will only consider the non- 
analytical contributions here, since we are only interested 
in non-local effects. 

As in the high energy renormalization of the the- 
ory will also be of no concern for us, since we are also 
only interested in finding the leading finite non-analytical 
momentum contributions for the 1-loop diagrams in the 
low energy scale of the theory. The singular analytical 
momentum parts which are absorbed into coefficients of 
the higher derivative couplings, will have no part to play 
here, they will ultimately not be manifested in this en- 
ergy regime of the theory. 

1. Defining the Potential 

The S-matrix is defined as the scattering matrix be- 
tween incoming and outgoing particles. The invariant 
matrix element iM. originating from the diagrams is 

(fclfc2 . . .| iTlkAks) = i2T:)H^{kA + kg - Efcfinal) 

(35) 

here we have two incoming particles. If we Fourier trans- 
form the earlier mentioned non-analytic terms to real 
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space, we easily see how the non-analytic terms con- 
tribute to the long-ranged corrections 



(27r)3 



(27r)3 



iq-r 



1 

Atkt 
1 

27r2r2 
1 



27rr'^ 



(36) 
(37) 
(38) 



obviously these terms indeed do contribute to the long- 
range corrections. When we calculate the tree diagrams, 
we explicitly see that the non-analytic contribution of the 
type (|36|l . will correspond to the Coulomb and Newtonian 
part of the potentials and the higher power of ^ will 
generate the leading order and classical corrections to 
the Coulomb and Newtonian potentials. Explicitly the 
invariant matrix element will look like 



"26 )^ 



) 2 2 2 

>2K e q 




(39) 



where A,B,... correspond to the analytical, local and 
short-ranged interactions, these terms will only domi- 
nate in the high energy regime of the effective field the- 
ory, whereas ai, a2, ■ ■ ■ and /3i, /32, . ■ . correspond to the 
leading non-analytic, non-local, long-range contributions 
to the amplitude. Many diagrams will yield pure an- 
alytic contributions to the S-matrix, such diagrams will 
not be necessary in our calculations, we will only consider 
the non-analytic contributions from the 1-loop diagrams. 
The diagrams which will yield non-analytic contributions 
to the S-matrix amplitude are those containing two or 
more massless propagating particles. 

Relating the Born approximation to the scattering am- 
plitude in non-relativistic quantum mechanics, we get in 
terms of iT 

(fcifca ...\tT \kAkB) = -iV{q){2TT)6{E - E') (40) 

where q = p' — p and V^(q) is the non-relativistic po- 
tential transformed in momentum space. We should be 
careful when comparing with (|35|l . in (iM) factors of 
(2mi X 2TO2) arise due to relativistic normalization con- 
ventions, thus we divide with these to obtain the non- 
relativistic limit. Equating the two we deduce 

-iVi(i){2TT)SiE - E') ^ {2TTf 5^{kA + kB - Sfcfi„al)f«>l 



or rather 



V{ci) 



1 I r (fk 



2mi 2m2 J (27r)3 



(41) 



2Try 6^kA+kB~^kRn.iKM) 
(42) 



Momentum integration yields the non-relativistic poten- 
tial 



^(q) = - 



1 1 

2mi 2m2 



M 



or in coordinate space 

y(x) 



1 1 



d^k 



2mi 2to2 J (27r) 



(43) 



(44) 



In our calculations, Ai will only contain the non-analytic 
contributions of the amplitude of the scattering process 
to 1-loop order, and we will not compute the full ampli- 
tude of the S-matrix, only the long-range corrections will 
be of our interest. In order to obtain their contribution 
to the potential, only a subclass of scattering matrix di- 
agrams will be required. If we wanted to find the full 
total non-relativistic potential, we would merely have to 
include the remaining 1-loop diagrams. This type of cal- 
culation has e.g. been performed in |23| (who also have 
used the same definition of the potential as us) where 
the full amplitude is considered. Their choice of poten- 
tial included all 1-loop diagrams, hence they obtained a 
gauge invariant definition of the potential. This choice of 
the potential makes good physical sense since it is gauge 
invariant, but other choices are also possible. The most 
convenient choice could depend on the physical situa- 
tion at hand or how the total energy is defined. The 
gauge invariant choice is also equivalent to the sugges- 
tion in (iJl , where it is suggested that one should use the 
full set of diagrams constituting the scattering matrix, 
from which one can decide the non-relativistic potential 
from the total sum of the 1-loop diagrams. However, it is 
worthwhile to note that we consider all the non-analytic 
corrections to 1-loop order, thus if we had the full am- 
plitude to 1-loop order we would still need to extract the 
non-analytical parts! We will continue using this defini- 
tion of the potential. 



V. RESULTS FOR THE FEYNMAN DIAGRAMS 

A. Diagrams contributing to the non-analytic parts 
of the scattering matrix potential 

In this section we shall extract the non-analytical parts 
of a limited set of 1-loop diagrams needed for the 1- 
loop scattering matrix in the combined quantum theory 
of QED and general relativity (however, it is a practi- 
cally complete set of diagrams in terms of non-analytical 
contributions to the scattering matrix). We will explic- 
itly see that the non-analytic contributions indeed corre- 
spond to the long range corrections of the potential. This 
will become obvious when the amplitudes are Fourier 
transformed to produce the scattering potential, whence 
all the analytic pieces are disregarded. The resulting 
non-analytic piece of the scattering amplitude will then 
be used to construct the leading corrections to the non- 
relativistic gravitational potential. 
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B. Classical physics 

Here we will look at the tree diagrams. The fermion- 
fermion scattering process at tree level should of course 
reproduce the results of classical physics both for gravi- 
tational interactions as well as for electromagnetic inter- 
actions. 



get when coupled to gravity. This difference will play 
a much bigger role when more complicated diagrams are 
involved. Indeed, the next set of diagrams were perhaps 
the most challenging of them all, the box and crossed box 
diagrams. 

C. The IPI diagrams 



1. Tree diagrams 

Given in figure ^ we have depicted the scattering pro- 
cess, where the (incoming/outgoing) momenta for parti- 
cle one are (k/k') with the (mass/charge) being (mi/ei), 
and similarly for the second particle with the (incom- 
ing/outgoing) momenta {p,p') and (7712/62) being the 
(mass/charge). This is assigned for all the other diagrams 
as well. The formal expression for the digram depicted 
in figure ^a), the scattering process involving a photon 
exchange, is 



We will calculate all the relevant IPI diagrams neces- 
sary to find the long range corrections to the potentials. 
We will start with the box and crossed box diagrams and 
continue with the set of triangular diagrams. Lastly we 
will work out the circular loop diagram. 



1. The box and crossed box diagrams 

There are in all four distinct diagrams. Two box and 
two crossed box diagrams, these are depicted in figure |21 
We will not treat all of these diagrams here. We will only 





1(a) 



1(b) 



FIG. 1: The set of tree diagrams 



iMna) = u{p')[T'']uip)u{k')[T^]u(k) 

and a graviton exchange, figure ^b) 

zXl(b) = u{p')[T^''Hp)u{k')[T'-''Hk) 



(45) 



(46) 



yielding the well known classical results, namely the 
Coulomb 



and Newtonian potentials. 

Vm{r) = - 



6162 
Attt 



Gmim2 



(47) 



(48) 



It is worthwhile to note that already at this stage the 
level of difficulty is not obvious. There is virtually no 
problem in working out the Coulomb term for the in- 
teraction, technically and mathematically it is straight- 
forward. However, in comparison, the Newtonian term 
is much more sophisticated to work out. This is mainly 
due to the many 7-matrices involved (since we are work- 
ing with fermions) and the long vertex rules theories 




2(a) 



2(b) 





2(c) 2(d) 
FIG. 2: The set of box diagrams contributing to the potential 

look at one of these, the rest can be treated in the same 
manner. Explicitly diagram |21[a) is defined by 



dH r 



(27r)4 



xu{p') rP^ip- t,p')DF{p- 1)t\p,p~ t) u{p) 

X u{k') Irt'^ik + £, k')DF{k + i)T^{k, k + I)] u{k) 

(49) 

The methods and techniques to work these diagrams 
are identical to those shown in we will repeat them 
shortly here. 

The only difference lies in the fact that these diagrams 
require four different integrals that were not worked out 
previously. We have worked them out, and the coeffi- 
cients can be obtained by contacting us, they are too 
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tedious to be written down in the appendix. Other than 
these integrals, these diagrams simply had to be worked 
out even though they involved enormous amounts of cal- 
culations. The level of difficulty is much larger than in 



J 



the previous case, due to the reasons mentioned earlier. 
All the box diagrams have been calculated by using sym- 
bolic manipulation on a computer. These have partly 
been checked by hand. 



On the mass shell we will have the following type of identities 



1 



£-p- 



-.{{l-vf-ml-e) 



(50) 



so 



{£ ■ q)r 



(27r)4 £2{e + qY [{£ + fc)2 - ml] [{I ~ pf ~ ml 



Si 

2 



d*£ 



£" 



(27r)4 £-2{£ + qf [{£ + kf - m^ [{£ - - m\ 



2 



(51) 



since the terms with {£ -f q)^ and simply do not contribute with non-analytical results. 

A more important reduction of the integrals is with contraction of the sources momenta rather than the exchange 
momentum 



d''£ 



{£ ■ fc)r 



(27r)4^2(£ + g)2[(£+y!j)2 

or in similar manner 

dH {£ ■ pY" 



m^] [{£ — p)2 — m: 



{2ny £'^{£ + q)^[{£ + ky - ml][ie - pY - ml 



d*£ 



£" 



(27r)4 £2{£ + qy[{£^p)2 - ml 



^ 2 P 



(52) 



dH 



£" 



(27r)4 £2{£ + q)^[{£ + ky - ml 



= -7;Ik (53) 



where the subscripts k and p on the /'s are written to indicate that the propagators left in the integrals are either 
from the particle with momentum k or p. Thus a loop momentum contracted with a source momentum simplifies our 
integrals considerably. 



The kinematics are (on the mass shell) 



2. The triangular diagrams 



1 I y 

k-q = p -1= Y 

q^ 

k ■ k' = m? 



k' ■ q = p ■ q 



p-p 



2 



Diagrammatically it is given in figure |21 these are 

(54) the only triangular diagrams that contribute with non- 
analytic terms. We will again only consider one instance 

(55) of these diagrams, namely the first fig. Ola). Formally it 
is written down as follows 



The potential contribution from these diagrams are found 
to be 



y{r)2(a)+2(b)+2{c)+2(d) 



3 eie2(mi 

4 Trr 



m2)G 



118Geie2 



(56) 



These diagrams yield both a classical contribution - the 
-pi - and a quantum correction - the ~ ^ . In the scalar 
QED calculations the box diagrams only generated quan- 
tum corrections, it is interesting to see that the Feynman 
diagrams do not necessarily generate identical results for 
the individual diagrams. This will become more clear in 
the other diagrams. 



iM 



3(a) 



d^£ 



u{p') 



u{p) 



xu{k') tI"'{-£- k,k')DF{-£- k)T'^{k,-£- k) u{k) 



mi 

£^ 



(57) 



Upon applying contractions and insertion of the relevant 
integrals where after Fourier transformations are per- 
formed, we end up with the potential contribution 

V3(a)+3(fc)+3(c)+3(d)(0 ^ ^^^^ 

which is also different from the result obtained in the @ 
for the triangular diagrams. 
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3(a) 



3(b) 





3(c) 



3(d) 



FIG. 3: The set of triangular diagrams contributing non- 
analytically to the potential. 



3. The circular diagram 
The circular diagram is depicted in figure 0] Formally 




4(a) 

FIG. 4: The circular diagram with non-analytic contributions, 
it can be written as 



(a) 



(27r)^ 



-Mp) 



fiupa 



u{p) 
u(k') 



-A"'(7) 



u{k) (59) 



When doing all the contractions and rearranging the 7- 
matrices one arrives at the result that the contribution to 
the potential from the circular loop diagram is precisely 
equal to nill. 



(60) 



Which indeed is very different from the result obtained 
in 0! 



D. The vertex correction diagrams 

There are several sets of IPR diagrams. All these are 
presented in this section. 



1. The IPR diagram 

The first set of these IPR diagrams is given in figure|Sl 
These diagrams are the only ones corresponding to the 





5(a) 



5(b) 



FIG. 5: The first set of IPR diagrams contributing non- 
analytically to the potential 

gravitational vertex corrections. Again we will only con- 
sider one instance of these diagrams. The matrix element 
originating from figure |SJa) diagram is 



{p') T^^DFip-i)T'' 



g2 



1^ 



u{k')[T^"'{k,k')]uik) (61) 



which yields the following potential when all the dia- 
grams are summed 



^5(a)+5(b)(0 = G 



(62) 



This checks with where it has also been calculated 

Of the next set of the IPR diagrams, depicted in fig- 
ure we will also only consider the first one. This is the 
first of the set of diagrams corresponding to the photonic 
vertex corrections. Formally diagram|^a) is given by 



iM 



6(a) 



(2^ 



{p') T'''DF{p-e)T^''{p,P~i) U{p) 



giving the potential 



[TP"^''-\q,l + q)] 



u{k')[T'']u{k) (63) 



6(a)-|-6(fc)+6(c)+6(rf) 



(r) 



36*6162 



(64) 
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6(a) 



6(b) 





2. The vacuum polarization diagram 

The diagram is depicted in figure IS] The formal ex- 
pression for this diagram is 



u{p')[t'']u{p) 



X T 



(Tp{5o 



-iri<i>e 



u{k')[T't']u{k) 



(67) 



This is the only instance in these calculations that 7- 



6(c) 



6(d) 



FIG. 6: The second set of IPR diagrams contributing to the 
potential 




7(a) 





7(b) 





FIG. 8: The vacuum polarization diagram contribution to the 
non-relativistic potential. 

matrices are not explicitly involved. Simple index con- 
tractions are done and one obtains the potential contri- 
bution after going to the non relativistic limit 



Geie2 



which is identical with the bosonic version. 



(68) 



7(c) 



7(d) 



FIG. 7: The third set of IPR diagrams contributing to the 
potential 



The second set of diagrams corresponding to the pho- 
tonic vertex corrections are given in figure \7\ Formally 
diagram [T^a) is given by 



(2^)4 

xuip')[T-P''^\p,p',e2)]u{p) 
X u{k')[T^]u{k) 

yielding the potential 



[e + qy 



y7{a)+7{b)+7{c)+7{d)ir) 



3Geie2 



(65) 



(66) 



VI. THE RESULTS FOR THE POTENTIAL 

When adding up all the separate contributions, we end 
up with 



V{r) 



Gmim2 ^ aeie2 
r r 

+ 2 (^"^261 + niie^ 



Gee eie2{mi + m2)aG 



4 Gah / 2^7*2 
3 TTc^r'^ V ^ mi 



2™l\ 

62 — 

7712/ 



15 



^ Gaheie2 



(69) 

where we have included the appropriate physical factors 
of c and we have further rescaled everything in terms 
of (5 = lastly (61,62) are the normalized charges in 
units of elementary charge. The result is divided into 
three separate parts, the first two terms represent the 
Newtonian and Coulomb potentials, the next two terms 
represent the classical post-Newtonian corrections to the 
potential, which also can be found by pure classical treat- 
ment of general relativity with the inclusion of charged 
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matter sources It is interesting to see that loop cal- 
culations also reproduce classical results, and not only 
quantum corrections. Finally, the last two terms are the 
leading 1-loop quantum corrections. We have in a greater 
extent reproduced the results of except for the last 
quantum correction where we get the factor (—15^) in- 
stead of (—8) as in Q. 

It is interesting to note how the different terms come 
about in the corrections to the potential. The classi- 
cal terms, the Newtonian and Coulomb terms, of course 
originate from the tree diagrams alone. The first post- 
Newtonian term and the first quantum correction no- 
tably originate from one type of diagrams alone, namely 
the gravitational vertex correction, see figure These 
match |Q] exactly. Now the second post-Newtonian cor- 
rection originates from all the box diagrams and nothing 
else, which indeed also is in full agreement with . How- 
ever, the second quantum correction is a sum of partly 
the box contributions and the rest of the diagrams. This 
correction does not match the quantum correction in 
ref. The rest of the diagrams have been done by 
hand, and all the diagrams done in this paper have been 
checked by symbolic manipulation on a computer. 



VII. DISCUSSION 

We have examined the leading order quantum cor- 
rections to gravitational coupling of a spin-i massive 
charged particle. Explicitly we have extracted the non- 
analytic terms from the diagrams, which exactly mani- 
fested themselves as corrections to the long range forces, 
this was realized when we Fourier transformed these 
terms into coordinate space. These terms originated from 
the propagation of the massless particles, here the pho- 
tons and gravitons. We have obtained similar results for 
most of the contributions to the corrections of the poten- 
tials, when compared with ^6j. Only one of the leading 
quantum corrections differs from the bosonic calculation. 

One could in similar manner do a QED - pure gravity 
scattering calculation as in However, more obstacles 
would be needed to overcome. First of all, one would have 
to find many new vertex rules involving the antisymmet- 
ric fields. Moreover we would have to derive second order 
Lagrangian densities both in terms of the symmetric and 
antisymmetric fields. Indeed, in this direction there lies 
a considerably interesting project ahead. 



APPENDIX A: FEYNMAN RULES 
1. Propagators 

The relevant propagators are presented in this section. 



a. Photon propagator 

The photon propagator is no stranger in QFT. In Feyn- 
man gauge the propagator becomes 

, _ -ir]a(3 



b. Graviton propagator 

The graviton propagator is perhaps a stranger. How- 
ever it has been worked out several places. In the Har- 
monic gauge we get the following for the graviton prop- 
agator 



q q-^ + ie 

with the projection operator 



+ ie 



c. Fermion propagator 

The fermion propagator can be found many places in 
literature, it is very well known 



k? — rn? 



2. Vertices 

The vertices are presented here. They are all derived 
in . For all vertices the rules of momentum conserva- 
tion has been applied. 



a. 1 -photon- 2- fermion vertex 
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The l-photon-2-fermion vertex can also be looked up 
in literature it is worked out to be 
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with 



6. l-graviton-2-fermion vertex 
The l-graviton-2-fermioii vertex is found to be 




where 



1 



c. l-photon-l-graviton-2-fermion vertex 

The l-photon-l-graviton-2-fermion vertex is not 
known from hterature, however it is found to be 



al3,k 



X 



{p,p') 



where 



1/3(7) 



{P,P') = ^7a(277^^"^-ry'^^"ry«") 



d. l-graviton-2-photon vertex 

We have derived the foUowing for the l-graviton-2- 
photon vertex 




o/3(7<5) 



{P,P') 



7,P 



where 



-paPifS) defined as above and {} do not involve any 
symmetrization factors, they merely exchange the in- 
dices. 



7,9 



APPENDIX B: TABLE OF RELEVANT INTEGRALS 



The following integrals are needed, note that in these integrals only the lowest order non-analytical terms are 
presented 

f dH 1 « r nrl 



together with 



327r2 
i 



(27r)4 P{l + qY 
(27r)4 £2(^ + ^)2 

(2^)4 £2(^ + ^)2 = ^ 



327r2 
i 



%qu[ - -L] ~q r]f. 



I = 



1 



{2tt)^ e'^{£ + q)^{{e + k)^ -m^) 327r2m2 
_d^ ijj, 

(27r)4 £2(^ + g)2((£+fc)2 „^2) 



[-L-S] 



327r2m2 



f f 



327r2m2 



q^Lql' [ ^ ^ ^ g*^ 



' 2 8 TO^ 



3 g2 



■1 , 1 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 



(B6) 
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iff 



(27r)4 P{£ + qf{{f + kf - m^) 
5 

16^ 



327r2m2 



1 



QiiQuqa [ L + —S ) + kf^k^ka y - g 

'(^Q/iki^ka H~ Qi/kf^ka + Qakfj,ki/^ I — + —— jS* 



3 TO"^ 16 m 
3 2to2 



+ {qf^q„ka + qt^Qaki, + qi^qakf, 

/I 1 



32 rn^ 



+ 



(B7) 



where 5* = ^ and L = ln(— g^). The ehipses denote higher order non-analytical contributions as well as the 

neglected analytical terms. Please note that there seems to be a typo in 6], in I^^a the factor after {k^k^ka) is 
lacking an L. Other than this typo all the integrals have been checked and are agreed upon. 
The following integrals are needed to do the box diagrams. 



K 
K' 







1 






(27r)4 £2(^^ 


hg)2((£- 


Ffc)2 -m2)((f - 


^p)2- 


- ml) 


d'^f 




1 






(27r)4 f^{e^ 


h<7)2((£- 


hfc)2 -m2)((£- 




- ml) 


d-^e 














-A:)2-m2][(^- 


p)2- 


21 ~ 


d^i 










(27r)4 




-fc)2-m2][(£-F 






d^t 










{2n)^ f^{l^ 


^qme^ 


-A:)2-m2][(^- 


p)2- 


to2] 



- 1 



3mim2 / 
W 



L 



3mim2 



L 



(B8) 
(B9) 
(BIO) 
(Bll) 



q^q^a + k^k^b + p^p^c + [qt^k^ + qyk^)d + [q^p^ + qi,p^)e + {Pf^k^ + Pukf^)f + i']^,„q'^g (B12) 



(27r)4 £2(^ + g)2[(£ ^. /,)2 _ ^2][(£ _^p,)2 _ ^2 



g^g.a' + fc^fc.fe' + p'^pIcJ'^ + {q^k, + q,k^)d' + {q^p^ + q,p'^)e' + (p'^k, + p',k^)}' + r/^,q'g' (B13) 
Here we have defined w = {k ■ p) — mim2 and = (fc • p') — miTO2- From these we can deduce an important relation 

_ 2 

that becomes vital during the calculations, namely W — w = k ■ (p' — p) — {k ■ q) — The w and W are displayed 

here only for the K and K' integrals, for the rest of the integrals we have used W = {k ■ p) — mim2 — see j23| for 
derivations. The coefficients to these are long and tedious to write down properly, however, if required, they can be 
obtained from us by contacting us. 

For the above integrals the following constraints for the non-analytical terms can be verified directly on the mass- 



shell. We will use hi • 



~ £2(f+q)2[(£+fe)2_^2] for k and p/ ~ j2(j+^yTpT^)2T^ 
no particular choice is needed for contractions with q 



for p, p>I ■ 



for p' and 



K n^^" — K' n'^'' — T n^^^ — T n^^" — J n^^" — Cl 



2 2 2 

^'^K K a'" — - — K K' n^^ — - — K' K' n" 



-K' 



^>'^ = ^fc^., Kp"'^\uI, K'^,k^^\p,h, KX^^p,I 
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J a ^ T T V ^ T T ^ ^ T J V ^ J 7M ^ J 



7/^ k°' — — T , T h'^ — - J , T ht^ — — T 

There seems to be a typo in 6] the metric in Ip.ua'H^'^ is written as I^ua'r]"'^ ■ 
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This is in analogy with the afhne connection relation F^^ 

that followed from DaQ^v ~ gtiv,a = 

The q is attached to distinguish the charges from the 

vierbeins. 



